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Abstract 

Bandyopadhyay and Das [Phys. Plasmas, 9, 465-473, 2002] have derived a nonlinear macroscopic 
evolution equation for ion acoustic wave in a magnetized plasma consisting of warm adiabatic ions 
and non-thermal electrons including the effect of Landau damping. In that paper they have also 
derived the corresponding nonlinear evolution equation when coefficient of the nonlinear term of 
the above mentioned macroscopic evolution equation vanishes, the nonlinear behaviour of the ion 
acoustic wave is described by a modified macroscopic evolution equation. But they have not 
considered the case when the coefficient is very near to zero. This is the case we consider in this 
paper and we derive the corresponding evolution equation including the effect of Landau damping. 
Finally, a solitary wave solution of this macroscopic evolution is obtained, whose amplitude is found 
to decay slowly with time. 
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I. INTRODUCTION 


Bandyopadhyay and Das derived a macroscopic evolution equation to investigate the 
nonlinear behaviour of the ion acoustic waves in a magnetized plasma consisting of warm adi¬ 
abatic ions and non-thermal electrons including the effect of Landau damping. This equation 
is a Korteweg-de Vries-Zakharov-Kuznetsov (KdV-ZK) equation except for an extra term 
that accounts for the effect of Landau damping. Bandyopadhyay and Das [l] reported that 
this macroscopic evolution equation admits solitary wave solution propagating obliquely to 
the external uniform static magnetic field and having a .sech 2 -profile. But the amplitude of 
solitary wave does not remain constant; it varies slowly with time r as (l + r/r')~ 2 , where t' 
is a constant depending on the initial amplitude of the solitary wave, the angle between the 
direction of propagation of the solitary wave and the external uniform static magnetic field 
and the parameters involved in the system. This evolution equation, which we are discussing 
about, losses its validity when the coefficient of the nonlinear term of the macroscopic evolu¬ 
tion equation vanishes and this vanishes along a particular curve in the /3<r-parametric plane 
as shown in Fig{T] where /3 is the nonthermal parameter associated with the nonthermal dis¬ 
tribution of electrons and cr is the ratio of the average temperature of ions to that of electrons. 
In this situation, in the same paper, they have derived a modified macroscopic evolution 
equation when the coefficient of the nonlinear term of the macroscopic evolution equation 
vanishes. This equation is a modified Korteweg-de Vries-Zakharov-Kuznetsov (MKdV-ZK) 
equation except for an extra term that accounts for the effect of Landau damping. Bandy¬ 
opadhyay and Das [l] reported that this modified macroscopic evolution equation admits 
solitary wave solution propagating obliquely to the external uniform static magnetic field 
and having a sech-prohle. But the amplitude of solitary wave does not remain constant; it 
varies slowly with time r as (1 + r/r') -1 , where r' is a constant depending on the initial 
amplitude of the solitary wave, the angle between the direction of propagation of the soli¬ 
tary wave and the external uniform static magnetic field and the parameters involved in the 
system. But again this modified macroscopic evolution equation is unable to describe the 
nonlinear behaviour of the ion acoustic waves including the effect of Landau damping if the 
coefficient of the nonlinear term of the macroscopic evolution equation approaches to zero, 
but not exactly equal to zero. In such situation, i.e., when the coefficient of the nonlinear 
term of the macroscopic evolution equation approaches to zero, but not exactly equal to 
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zero, Das et al. [2| derived a combined MKdV-KdV-ZK equation to describe the nonlinear 
behaviour of the ion acoustic wave when the Landau damping effect has not been taken 
into account. Following Ott and Sudan j^|, in the present paper, we derive a macroscopic 
evolution equation to study the nonlinear behaviour of the ion acoustic waves including the 
effect of Landau damping. This equation is a combined MKdV-KdV-ZK equation except for 
an extra term that accounts for the effect of Landau damping. The solitary wave solution 
of this further modified macroscopic evolution has been obtained, ft is found that due to 
inclusion of the effect of Landau damping the amplitude of the alternative solitary wave 
solution of this equation is a slowly varying function of time. 


In the investigations made by Das et al. 


s, 


the Landau damping effect has not been 


taken into account. In the present paper, we include this effect on the problem considered 


in Das et al. [2|. Starting from the same governing equations but replacing the expression 
for the number density of non-thermal electrons by the Vlasov-Boltzmann equation for 
electrons, an appropriate macroscopic evolution equation corresponding to the combined 
modified Korteweg-de Vries-Zakharov-Kuznetsov (combined MKdV-KdV-ZK) equation of 
Das et al. Q] is derived, which describes the long-time evolution of weakly nonlinear long 
wave-length ion acoustic waves in a magnetized plasma consisting of warm adiabatic ions 
and non-thermal electrons including the effect of Landau damping. 

The physics of nonlinear Landau damping is of interest for two major reasons. First, it is 
a fundamental and distinctive plasma phenomenon that links collective and single-particle 
behaviour. Second, the derivation of reduced fluid models that incorporate accurately such 
kinetic effects, is of great importance for plasma transport studies. For instance, some 
authors have proposed a fc-dependent dissipation term, which correctly reproduces linear 
Landau damping within the framework of fluid models 4]. However, the long time behaviour 
of Landau damping is intrinsically nonlinear, and, in order to assess the validity of the above 
models, it is important to understand whether the damping will continue indefinitely, or will 
eventually be stopped by the nonlinearity. 

The research works on solitary waves in plasmas have been done under various physical 


colons such as plasmas including muh.-spemes ions & negative ,ons 0, and dust pan 
tides [7j]. In many cases, the Korteweg-de Vries (KdV) equation is used to describe basic 
characters of the wave. Detailed properties of the solitary waves observed in experiments 
in plasmas are, however, slightly different from those predicted by the equation. Using 
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Q-machine plasmas, Karpman et al. jg] have observed oscillations in the tail of solitary 
waves, which are caused by resonant particles and have shown that the tail changes its 
shape depending on the strength of Landau damping. 

This paper is an extension of the work of Bandyopadhyay and Das jl| , where we derive a 
further modified macroscopic evolution equation which describe the non-linear behaviour of 
ion-acoustic waves in fully ionized collisionless plasma consisting of warm adiabatic ions and 
non-thermal electrons having vortex-like velocity distribution, immersed in a uniform static 
magnetic field directed along z-axis including the effect of Landau damping. This equation is 
true only for the case when the coefficient of the nonlinear term of the macroscopic evolution 
equation derived by Bandyopadhyay and Das [lj] approaches to zero and not exactly equal 
to zero. With the help of multiple time scale analysis of Ott and Sudan Q, we find a solitary 
wave solution of this equation. From the solution, we can conclude that the amplitude of 
the solitary wave slowly decreases with time. 

This paper is organized as follows. The basic equation have been given in Section [TT1 
The macroscopic evolution equations are given in Section m in which the derivation of 
MKdV-KdV-ZK like macroscopic evolution equation is given in subsection IIII Cl Solitary 
wave solutions of the combined MKdV-KdV-ZK like macroscopic evolution equation are 
investigated in Section ITVl Finally, we have concluded our findings in Section fVl 


II. BASIC EQUATION 


The following are the governing equations describing the non-linear behaviour of ion- 
acoustic waves in fully ionized collisionless plasma consisting of warm adiabatic ions and 
non-thermal electrons having vortex-like velocity distribution, immersed in a uniform static 
magnetic field directed along z-axis. Here it is assumed that the plasma beta i.e., the ratio 
of particle pressure to the magnetic pressure is very small and the characteristic frequency 
is much smaller than ion cyclotron frequency (Cairns et al. 9], Mamun 0 )- 


dn 

dt 


+ V.(nu) 


0, 


( 1 ) 


<9u 

~dt 


(u.V)u = —'Vy? + w c (u x z) 


er 
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-Vp, 


( 2 ) 
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X7 2 ip = n e — n, 


( 3 ) 


p = n 1 . 


(4) 


where 


n P = 


fdv n , 


(5) 


and the velocity distribution function of electrons / must satisfy the Vlasov- Boltzmann 
equation 


frn^df_ df_ &pdf_ = 

m,- dt V " dz dz dv n 


( 6 ) 


Here n,n e ,u,p,ip,(x,y,z) and t are respectively the ion number density, electron number 
density, ion fluid velocity, ion pressure, electrostatic potential, spatial variables and time, 
and they have been normalized respectively by no (unperturbed ion number density), no, 
c s ( = \J ) (ion-acoustic speed), n 0 K B T h A D ( = yj ) (Debye length) and 
a;" 1 (ion plasma period), where a = is the ion cyclotron frequency normalized by 


oj p \ = ■ /^znioe 


and y( = |) is the ratio of two specihc heats. Here K B is the Boltzmann 
constant; T e ,Tj are respectively the electron and ion temperatures; m is the mass of an ion 
and e is the electronic charge v„ is the velocity of electrons in phase space normalized to 


v, = 


■\J K .^ e • In (SJ), the adiabatic law has been taken on the basis of the assumption that 
the effect of viscosity, thermal conductivity and the energy transfer due to collision can be 
neglected. 

Since the electrons are assumed to be nonthermally distributed, the electron velocity 


distribution function can be taken as (Cairns et al. 0 ) 

1 


/oK) = 


(l + aiu,f)e ( 2 ^) 


(7) 


\/27r(l + 3cq) 

To discuss the effect of Landau damping on ion-acoustic solitary waves, we follow the 
method of Ott and Sudan [3] and following them, we replace by •> /IF £ where £ is a 
small parameter. The equation (|6]) then assumes the following form 


df df dip df 

a *dt +v "dz + dza^r °’ 


( 8 ) 


where a 2 = 


5 













(9) 


Again using Eq. (j3J), the equation (J2]) becomes 

— + (u.V)u = — \7ip + cu c (u x z) — crnsVn, 
ot 

and the equation (J3]), which is the Poisson equation becomes 


OO 


+ n 


fdv„ = 0, 


( 10 ) 


—OO 


Therefore Eqs. (JT]), (J9]), (TTOl) and (jHJ) are our governing equations. 

III. EVOLUTION EQUATIONS 
A. Macroscopic evolution equation 

Before deriving the nonlinear evolution equation for ion-acoustic wave in a magnetized 
collisionless plasma consisting of warm adiabatic ions and non-thermal electrons including 
the effect of Landau damping for a particular case not considered in the paper of Bandy- 
opadhyay and Das [l|, we give below in short a summary of the results obtained in that 
paper. The macroscopic evolution equation obtained is the following: 



( 11 ) 


where 



( 12 ) 



(13) 



( 14 ) 
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(15) 


E = 


1/ 


4\/27r 


(4 — 3/3), 


and the constant V is given by 


(i-/3)h/ 2 


-o = 1. 


(16) 


The equation (fTIj) is a KdV-ZK equation except for an extra term (last term of the left hand 
side of (TTT1) ) that accounts for the effect of Landau damping. The solitary wave solution of 
the equation (TITj) has been obtained in that paper of Bandyopadhyay and Das [1]- They 
have found that the solitary wave solution of the equation (1111) has the same seek 2 - profile 
as in the case of KdV-ZK equation. But, here the amplitude as well as the width of the 
solitary wave varies slowly with time. In particular, the amplitude (a) of the solitary wave 
solution of the equation (TlTT) is given by the following equation. 


a — cio( IT 


t 


-2 


(17) 


where a 0 is the value of a at r = 0 and r' is given by the following equation 
1 


r = 


, A 2 E 2 a 2 B'a 0 
4 \l 6(cos 2 5 + D sin 2 S) 


p r r S rfT (srfz,) dz ' dz 

J —oo J — oo 


dZ' 


Z-Z' 


-1 


(18) 


Using (pup , the expression for B' can be simplified as 

y 20(1 -Pf 


B' = 


9 


\o-o. 


/s)> 


(19) 


where 


9{1-3(1-OT 
40(1 — f3) 3 


( 20 ) 


From this expression of B', it is easy to see that the coefficient of the non-linear term 
in (1TTP vanishes along a particular curve Fig{T]in the /3o plane, and consequently, it is not 
possible to discuss the nonlinear behaviour of ion acoustic wave including the effect of Landau 
damping with the help of Eq. (fill) . In this situation, i.e., when B' = 0, Bandyopadhyay 
and Das |l] have also derived a modified macroscopic evolution equation. 
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B. Modified macroscopic evolution equation 


For this case, i.e., when B' = 0, giving appropriate stretching of space coordinates and 
time, and appropriate perturbation expansions of the dependent variables Bandyopadhyay 
and Das [l] in the same paper have derived the following modified macroscopic evolution 
equation for ion acoustic waves in a fully ionized collisionless plasma consisting of warm ions 
and non-thermal electrons immersed in a uniform static magnetic field directed along the 


z-axis: 


El + , l AD l (EEL , ay il f 

Sr + ^ ‘ SC + 2 A SC 3 + 2 AD dc\ sc 2 + djf J 


— AE0I2P 


dyV d(' 


= 0. 


SC' c - c - (21) 

Here A, D and E are same as given by the equations (fl2l) . (TT4l) and (1T5l) respectively and 
B" is given by the following equation: 


B" = 


V 2 - a 


:(1 + 3/3) 


-4 


\ 27 2 


V 2 - -a 


-1 


3V 2 - -a 
9 


( 22 ) 


and the constant V is determined from equation (1151) . 

The equation (l2Tj) is a MKdV-ZK equation except for an extra term (last term of the 
left hand side of (12ID ) that accounts for the effect of Landau damping. The solitary wave 
solution of the equation (ITT]) has been investigated by Bandyopadhyay and Das [l] in the 
same paper. They have found that the solitary wave solution of the equation f|2Tjl has the 
same sec/i-profile as in the case of MKdV-ZK equation. But, here the amplitude as well as 
the width of the solitary wave varies slowly with time. I 11 particular, the amplitude (a) of 
the solitary wave solution of the equation (j2TT) is given by the following equation. 


a — Uq ( 1 H —1 


-1 


where a 0 is the value of a at r = 0 and r' is given by the following equation 


r = 


A 2 E 2 a 2 2 a 2 0 B" 


1 c 

2 V 3(cos 2 5 + D sin 2 5) 


P 


sechZ 


d(secliZ') dZ'dZ 


1 -1 


' —OO J —OO 


(23) 


(24) 


8Z< Z-Z' 

But both the evolution equations (HIT) and EJ are unable to describe the nonlinear 
behaviour of the ion acoustic wave along with the effect of Landau damping in the neigh¬ 
bourhood of the curve o = ap in the /3a - parametric plane along which B' — 0 (Fig{l]). 
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This is the situation we are considering here. We have derive in this case, a further modified 
macroscopic evolution equation which describes the nonlinear behaviour of the ion acoustic 
wave in the neighbourhood of the curve a = Gp in the /3a - parametric plane along which 
B' = 0. 


C. Further Modified Macroscopic evolution equation 


To discuss the nonlinear behaviour of the ion acoustic wave in the neighbourhood of the 
curve in the (3a parametric plane along which B' = 0, we assume B' cs 0( e ^) (Nejoh [11]), 
and we take the following stretching of space coordinates and time. 


£ = £ 2 x, rj = £ 2 y, ( = e 2 (z — Vt), r = e 2 t, 


(25) 


where e is a small parameter measuring the weakness of the dispersion and V is a constant. 

With the stretching given by (J25]), the equations (JT]), (J9]), (TTOl) and (J 8 ]) respectively assume 
the following form: 


i dn 3 dn i 
-£ 2 V— + £2__ + £ 2 V ? .(mi) = 0, 
o( or 


(26) 


i Tr <9u 3<9u i. _ \ _ 5 i _i 

-£ 2 V—+£ 2 — —h £2 (u.V|)u = — £2 + o; c (u x z) — -a£ 2 n 3 V ? n, (27) 

U C, C/'T O 


f c 

eS7\ip + n- 

J — C 


fdv M = 0 , 


(28) 


Here 


T r s df s df i df i dip df 

~ VOL<2£ 2 ~77Z + Oq£ 2 77-1" £2V "~x7 + £2 ~^7 "a— — 0- 

df or df df dv„ 


_ ^ d ^ d ^ d 

= t. jp + e,Q- + e.-tp U = hi. v, to). 


(29) 


(30) 


Next we use the following perturbation expansions of the dependent variables to make a 
balance between nonlinear and dispersive terms. 


(' n,(p,w,f) = ( 1 , 1 , 0 , 0 , / 0 ) + J2Zi £2 ( n(i) i<P {i \ w{i) J (i) ), 
(u,v) = 


( 31 ) 
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Substituting (l3Tf into the equations (l26]) - (j29j) and then equating coefficient of different 
powers of e on both sides, we get a sequence of equations. From the lowest order equations 
obtained from 0261) - (1281) , which are at the order e, we get the following equations. 


77,(d = ( V 2 — -Irr 1 in W 


-1 

3 Uj 


W W = v(v 2 - for^ 


u (D = _Zl(y2_| 0 . 


dipW 
dr) ' 


v w = ¥i[v 2 -la 


« (1) = fZc / (1)rf br 


-1 


dip h) 

d( ■ 


(32) 


From the Vlasov equation ()29l) at the lowest order, i.e., at the order e 1 / 2 , we get the following 
equation 


df w a/,y) 

<9C (9t’„ <9<C 


0. 


(33) 


As this equation does not have a unique solutions, we include an extra higher order term 
= ££ 2 ol 2 ^ 7 - originated from the Vlasov equation at the order e 2 . Let us write the 
equation (1331) as follows. 


2 df t 
£ Oi 2 


(1) 


dr 


+ hi, 


df t 


(i) 


+ 8fo V 11 = „ 


d( dv„ d( 


(34) 


Then /« can be obtained as unique solution of this equation by imposing the natural 
relation of the form 


/ (1) = lim/ £ (1) - 


(35) 


Assuming r dependence of /i 1 and to be of the form expiiur ), the equation (j35j) can 
be written as 


iue 2 a 2 f^ + 


SfP 

S( 


+ 


dfo d<p {1) 
dv n d( 


(36) 


Now taking Fourier transform of this equation with respect to the variable ( according to 
the definition, 

g = g(k) = —!= [ g(()exp(-ik()d (, (37) 

V J —00 
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we get 


iuj£ 2 aofi 1 ' 1 + ikv,J ^ = 0. 

This equation gives the following expression for 


(38) 


fP = 


/ r Q/o 

/V' r> 

OV II 




(i) 


(39) 


kv„ + a 2 cu£ 2 

Now whenever the factor 1 /(kv„ + a 2 (jj£ 2 ) comes under an integration over v u along the real 
axis, the general prescription is to replace this integration according to Landau, along a 
contour in the complex u M -plane known as Landau contour. This is equivalent to replacing 
the factor l/(kv u + a 2 uj£ 2 ) by the following 

1 1 


= P- 


+ in 5(kv n + a 2 u£ 2 ). 


(40) 


kv I, + a 2 u>£ 2 kv„ + a 2 u£ 2 

Substituting this relation into the equation (1301) and then proceeding to the limit £ —y 0+, 
we get according to (1351) the following expression for / (1 b 

/ (1) = " 2 ^^ (1) “ 2n ^kv u S(kv u ). (41) 

Due to the relation x5(x) = 0, this equation assumes the following form: 


Taking Fourier inverse transform of 


dv 2 

, we get 


/<» = -2^. 


(42) 


(43) 


Substituting (1431) in the last equation of (f32l) and then performing the integration we get 

4a i 


n 


(b = (1 - p)^ 1 ) with P = 


(44) 


This equation along with the first equation of 
determine the constant V 


1 T 3cq 

gives the following dispersion relation to 


(l-P)[V 2 ~p\ =1. 


(45) 


This equation is same as the equation (ITBll as well as the equation (20) of Das et. al. 

In the next order, i.e., at the order £ 3 / 2 , solving the ion continuity equation and the parallel 
component (i.e., the component parallel to the ambient magnetic held, i.e., z-component or 


11 










(■-component) of ion fluid equation of motion for n^ and w^ to express them in terms of 
and ip( 2 \ we get the following equations: 


n (2) _ 


W& = 



3V 2 -Id 


<^ (2) + 


v2 -l< 


-tpW + \v- 


V 2 + f c 


V 2 -lc 


M‘>] 2 , 




(46) 


From the perpendicular component (i.e., the component perpendicular to the ambient mag¬ 
netic field, i.e., the components along x-axis and y-axis) of the ion fluid equation of motion 
at the order £ 3 / 2 , we get the following equation: 

du^ dv ® V 3 ( 2 5 \ _1 <9 /<9 2 </?W d 2 Lp < ^ , \ , . 

<9( drj cu 2 \ 3 / d( \ d ( 2 dp 2 )' 

From the Poisson equation at the order e, we get 


n 


( 2 ) _ 



0. 


(48) 


To find /l 2 ), we again consider the Vlasov equation at the order e 3 ^ 2 . The Vlasov equation 
at the order e 3 ! 2 is the following, in which as mentioned in the lowest order Vlasov equation, 
an extra time derivative term = £ 3 ^ 2 £ 2 ct 2 ^— has been included and f 1 - 2 ' 1 has been 

replaced by /J 2 . 


2 9f t 

£ 01 2 


( 2 ) 


dfP , dfodipW 


dr 


+ v "^ + -a£, — ~ 2v " 


a 2 h a 

d(v 2 ) 2 d( 




(49) 


Then /< 2 > can be obtained as unique solution of this equation by imposing the natural 
relation of the form 




lim fi 2 ( 

e->0+ j£ 


(50) 


Assuming r dependence of /J 2 ' and (p^ to be of the form exp(iujr), taking Fourier transform 
of this equation with respect to the variable (, using the causality condition (|40|) and finally 
proceeding to the limit e —> 0+, we get according to (15(1 the following expression for /^ 2 h 


/' 2) 




+ 2 


d 2 f 0 

d(v 2 ) 2> 


(51) 
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where 


d a = [0 (1) ] 2 , (52) 

and we have used the relations xS(x) = 0 and xP(^) = 1 to simplify the equation (T5T11 . 
Taking Fourier inverse transform of (|5TT) . we get 


f (2) = -2 




(53) 


dv 2 d(v 2 ) 2 

Substituting for /( 2 ) and given respectively by (T53T) and the first equation of ()46l) into 
the equation (1481) . we get the following equation after simplification 


- V 2 - -a 


-1 r 


(1-/3) (V 2 ^-a -1 


<p® + B'[<pM] 2 = 0. 


(54) 


Now the first term of the left hand side of the equation ()54l) is identically equal to zero due 
to the dispersion relation as given by the equation (1451) and as B' ca O^ 1 / 2 ), the second 
term of the left hand side of the equation (T54|) along with its sign has to be included in the 
next higher order Poisson equation, i.e., this term along with its sign must be included in 
the left hand side Poisson equation of order e 3//2 and consequently the equation (|54|) . i.e., the 
Poisson equation of order e 2 / 2 (= e) is identically satisfied. So, including the term +S'[(^d)] 2 
in the Poisson equation of order e 3 / 2 , we can write the Poisson equation at the order e 3 / 2 as 
follows 


d y 1 ) ay 1 ) d y 1 ) 


+ 


+ 


+ n 


(3) 


f (3) dv„ + B'[^ 1) ] 2 = 0. 


(55) 


<9£ 2 dif d( 2 

Now, in the next order, i.e., at the order e 2 , solving the ion continuity equation and the 
parallel component of ion fluid equation of motion for the variable to express it in terms 
of we get the following equation. 


l/ 4 


1 ty 3 ) 


sc + 2 T 2 -r 


-2 


yd) 

dr 


d 


CU 2 [ V 2 - fcx 


d( 


ay 1 ) ay 1 ) 


de 


dr] 2 


V 2 


-a 


-4 


10 , /2 3« 3l/2 -P 

—u — 6V +- 


V 2 


7<r 


[^ (1) ] 


2 y (1) 


3V2 ~I°, g 


v 2 -k 


d( 


[^ (1 V 2) L 


(56) 
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where we have used equations ([32]) . (1461) and (1471) to get this equation in this present form. 

Now our task is to find that determines n^ from the Poisson equation (l55j) at the 
order e 3//2 . To hnd /( 3 ) we consider the Vlasov equation at the order e 2 . The Vlasov equation 
at the order e 2 is the following, in which as in the lowest order case an extra higher order 
term e 4 a 2 ^^ - has been included and f (:y> has been replaced by /i 3 ' ) and where we have 
substituted the expressions for /d) and f^ 2 - given by equations (l43j) and (j53l) respectively. 


2 9f t 
e a 2 


(3) 


dr 


, dfr , dfpdcpW 

d( dv u d( 


ota d fo , , d2 fo 
= -2Va 2 — I x 2 + 4u„—— 

ov„ o(v t 


2\2 


y 2 - 4t’, 


d 3 fo 

"d{v 2 y 


r^2, (57) 


where 


x 2 


dip^ 


Z2 = [v m]^. _ 


(58) 


Therefore /^) is obtained from the unique solution of the equation (HTTP by the relation 


/ 


(3) 


lim /i 3) - 

£ —^0-K 


(59) 


As in the earlier cases, assuming r dependence of /j 3 ^ and <^ 3 ) to be of the form exp{iujT ), 
taking Fourier transform of this equation with respect to the variable (, using the causality 
condition (l40j) and finally proceeding to the limit e —> 0+, we get according to (150|) the 
following equation determining /®: 


ik 


/(»> + 


9 T/ d f ° 

2Va2 dv 2 
+4 n , o^ V2 


kPy-^—J + insgn(k)S(v„ 

A - 

4 n.^ 2- 


x 2 


Integrating 


d(v 2 y yz d(v 2 r 

over the entire range of n„, we get the following equation. 


(60) 


ik[n t 


(3) 

e 


(1 -/3)^ (3) ] = ~iVa 2 (i-3P)j£sgn(k)x2 + y2 + i(l + 3/?)i 2 , 


(61) 


where we set 



( 62 ) 
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Taking inverse Fourier transform of the above equation, we get 


dn t 


(3) 






V 


4v / 27t 


« 2 (4-3 P)P 


dcpW d(' 

d? C-C" 


(63) 


in which the convolution theorem has been used to find the inverse Fourier transform of 
sgn(k)x 2 - Now using the equations (|62l) and (1631) . we get the following equation. 


d_ 

5C 


/ (3) ch 


= ( 1 -«¥ + I ( ^) + 5 (1 + 3/ ’)[ vW ] > ¥ 
1 v^-wpr^ < 


4\/27r 


ac' C-C'’ 


(64) 


Substituting (j64|) into the equation obtained by differentiating the Poisson equation 
the order e 3 / 2 with respect to £, we get the following equation 


at 




d 


d( l ae 
a 


<9 2 (p^ 5 2 (pd) <9 2 <pd) 


ch? 2 


<9C 2 


+ 


4\/27r 


\/«: 2 (4 - 3/3)P 


<9<p« dc' 

ac ; C-C' 


— ((p (1) <p (2) ) = o. 

d( 


(65) 


Now substituting for given by (1561) into the equation fj65h . we get the following 
further modified macroscopic evolution equation, where the term +2B'ip^) being of 
higher order since B' = 0( e1 ^ 2 ) has been omitted. 


dipW , (1) dip (1) ,, r 1 . <9 3 <p (1) 


S( 

d ( 9 2 (ph) ^ 9 2 (pP) 
dr } 2 


2 AD 3(\ d£ 2 


+ 


+ —AEQ.2P 


d ( 3 

ac ; C-C 


= 0 . 


( 66 ) 


ffere A, B', D and P" are respectively given by the equations (fT2j) - (fT51) and (1221) and the 
constant 14 is determined by the equation (1T6)) . The equation (1661) is a combined MKdV- 
KdV-ZK equation except for an extra term (last term of the left hand side of (1661) ) that 
accounts for the effect of Landau damping. In the next section, we find the solitary wave 
solution of this further modified macroscopic evolution equation. 
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IV. SOLITARY WAVE SOLUTION OF THE FURTHER MODIFIED MACRO¬ 
SCOPIC EQUATION 


If we neglect the electron to ion mass ratio, i.e., if we set a 2 = 0, the equation 


reduce to a combined MKdV-KdV-ZK equation. The solitary wave solution of this combined 
MKdV-KdV-ZK equation has been studied in Das et al. [2]. In this paper, our aim is to 
find the solitary wave solution of the equation ()66|) . 

The solitary wave solution of the equation (l66j) with a 2 = 0 propagating at an angle 
5 with the external uniform static magnetic field is the following, which has already been 
obtained in section IV of Das et al. pj] , 


=<fo(Z) = o^. 


(67) 


where 


S = sech[2pZ ], 


( 68 ) 


f = S + A Vm, 


(69) 


A = ±1, 


(70) 


a = 


12p 2 (cos 2 5 + D sin 2 5) 
~B' 


(71) 


M 


12p 2 T?"(cos 2 5 + D sin 2 5) 

+ ga 


(72) 


Z = £ sin S + ( cos S — Ur. (73) 

For the existence of the solitary wave solution fl67]h it is necessary that the following condition 
is satisfied. 


L = MB' 2 = B' 2 + 125 V(cos 2 S + D sin 2 5) > 0. 
If the condition (1741) holds good, U is given by the equation 

U = 4p 2 a 3 , 


(74) 


(75) 
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where 


a 3 = -Acos<5(cos 2 5 + .Dsin 2 <5). 


(76) 


With the help of the equations (TTTTh (172lh d75l) and d76l) . we get the following expressions of 
p, U and M to express them in terms of a. 


v — — 
2 


aai 
6a 3 ’ 


U = —aai, 
o 

M — 1 a. 


a 2 
a 3 ' 


where 


(77) 

(78) 

(79) 


a 3 = AB'cosS , (80) 

a 2 = AB" cos <5. (81) 


Using ([77J), we can write the Eq. flfiTl) as 


<p (1) = ipo(Z) 


= a- 


seek 

sin 5 + ( cos 5 — |aair) 


seek 

sin <5 + C cos 5 — |aair) 

+ Ay/l + 


(82) 


Assuming that a to be a slowly varying function of time, following Ott and Sudan Q, we 
introduced the following space coordinate in a frame moving with the solitary wave. 


Z = 


—- f £ sin5 + (cosS - a± I adr ). 

6a 3 \ 6 


(83) 


It is important to note that if a is a constant, then Z = 2 pZ and consequently, 


<d (1) = <Po{Z) 

sechZ 

= a -—- 


sechZ + \VM 


seek 

(£ sin 6 + C cos 5 - |ai f Q r adr) 


seek 

\J sin <5 + ( cos 5 — / Q r adr) 

+A v /nw 


(84) 
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is the solitary wave solution of the combined MKdV-KdV-ZK equation propagating at an 
angle 6 to the external uniform static magnetic held. Now dropping “overline” on Z, we can 
write the equation (1841) as 


=MZ) 


sechZ 


= a 


sechZ + \VM 


= a- 


seek 

(£ sin 5 + C cos 5 — |«i f Q r adr) 


seek 

^/||^(£ sin 5 + ( cos 5 — |cq JJ adr) 

+ Ay/l + 


(85) 


where Z is given by the following equation: 

Z = 


aa i 
6a 3 


£ sin5 + (cos5 -aq 

6 


adr 


( 86 ) 

is a solitary 


Now our aim is to find the condition for which po> given by the equation 
wave solution of the further modified macroscopic equation (l66jh 

With the change of variable defined by the equation (1861) and assuming that yd 1 ) is a 
function of Z, r only, Eq. (166|) can be written as 


dp^ 

dr 


Z da\ (9yd 1 ) 


<9ydd 


+ ( — -cq P a + 7^—+ 2 pcq^ (1 ) + 2 pa 2 (p { ' l) ) 


d 3 p C 1 ) 


+8p 3 a 3 "d (1 + AEa 2 p cos 5P 

dZ 6 


2adrJ dZ ' " ±r dZ 
dp w dZ’ 


, dp^ 

dZ 


= 0. 


(87) 


_ ^ dZ Z-Z' 

To investigate the solution of Eq. (1871) . we follow Ott and Sudan [3|] and generalizing the 
multiple-time scale analysis with respect to a 2 , by setting 


pW(Z, r) = q<!» + a 2 qW + a \q^ + a\q ® + 


( 88 ) 


where each q^\j = 0,1, 2, 3,....) are the function of r = r 0 , Tj, r 2 . Here r 3 is given by 


Tj = a J 2 r,j = 0,1,2,3, 


(89) 


Substituting (1SS1) into (187| and then equating the coefficient of different power of a 2 on each 
side of Eq. fl8?j) . we get a sequence of equations. The zeroth and the first order equation of 
this sequence are respectively, given by the following equations. 


P 


' d_ 


Z da d 
2 a dr dZ 


M 


+ 4 « (0) 


0, 


(90) 
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p 


d Z da d 

d^ + 2 ad^dZ 


q (1) + 4? Lq{1) = p M( i {0) > 

dZ 


where 


i= 92 


+ - <Z (0 > + h^_il[ ? (0> r _ 1, 

dZ 2 a a 2 


(91) 


(92) 


C / 6a 3 -3 

P = !q a 3j 


af 


(93) 


Mg (0) = 


-7- 1 ----hAEpcosflP / -f- 

dr, 2a 5n dZ 1 J „ dZ' Z - Z' 


(94) 


sechZ 


Now it can be easily verified that g^ 0) = a ~ sechZ+x ^jj 
order equation if 


da 


= 0, 


is the soliton solution of the zeroth 


(95) 


which implies that a is independent of time, i.e., at the lowest order, the solitary wave solu¬ 
tion of the further modified macroscopic evolution equation is same as that of the combined 
MKdV-KdV-ZK equation. 

Using ()95lh Eq. (T9TT) can be written as 


p %r + ik Lq{ " 


(96) 


Now for the existence of a solution of the equation (l96]h its right hand must be perpendicular 
to the kernel of the operator adjoint to the operator L ; this kernel, which must tend to 
zero as \Z\ —» oo is ■ Thus we get the following consistency condition for the 

existence of a solution of the equation (J96J). 

f°° sechZ 


Mq {0) dZ = 0. 


(97) 


J -oo sechZ + \VM 

From equation (197)) . we get the following differential equation for the solitary wave amplitude 
a. 


da AEa 3 / 2 (B' + aB") cos 5 
dr\ ^/3P'(cos 2 5 + D sin 2 5) 


r»00 r oo 


P 


_r . r~7\ d r t / rv! M dZ'dZ 

<P(Z) — = 0, 


' —oo J —oo 


'Z-Z' 


where 


mz) = 


sechZ 


sechZ + \VM 


(98) 


(99) 
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Using the relation 77 = « 2 T, the equation 
form: 

da AEa 2 a 3 ^ 2 (B'+ aB") cos 5 


can be written in the following simplified 


POO /*oo 


p 


m ,^ mz'))dz'dz nnm 

dr ' ^3B'(cos^ + Dsin 2 i) ‘ •/-«, J-oo ( } SZ' Z - Z' ' 

Here it is important to note that M(— 1 + a^) appearing in 'i’(Z) is a function of a. So, 


it is not possible to find the exact analytical dependence of a on r. But we can solve the 
above equation by using the Taylor series expansion for the terms of the form i n 

powers of a. Keeping terms upto the order a 5//2 , we get the following differential equation 
for a from equation ( 1 1001 ) . 

da 


dr 


AEa 2 a 2 , 


A TP 

AEa 2 a 2 — 


B' 

3(cos 2 5 + D sin 2 5) 
B" 


2 5 n COs5 71 

cos S (72 + 73 ) 


AEa 2 a 2 


2 a/3H'(cos 2 5 + D sin 2 5) 

B" 

cos 5 71 = 0, 


a/3H'(cos 2 5 + D sin 2 5) 
where 71 , 72 , 73 are given by the following integrals. 


( 101 ) 


72 = PfZfZ^)Sl^(Z')}^, > 
Ts^nroo^iMfi, 

<E>i(Z) and <h 2 (Z) appearing in the above are given by 

sechZ 


( 102 ) 


$1 (Z) = 
$ 2 (Z) = 


sechZ + A’ 
sechZ 


(103) 

(104) 


(sechZ + A ) 2 

Now solving the above differential equation (11011) for a by the use of the initial condition, 
a = a 0 when r = 0 , we get the following equation for a: 


/i tan 


-1 


1 + n 2 y/aao 


y/a- y/a^ _ r 


/aao 


T 


where 


T = 


1 ATP I B ' 

2 a2 V 3(cos 2 5 A D sin 2 5) 


cos S 71 


-1 


(105) 


(106) 
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(107) 



_ A (72 + 73 ) 
2 


7i 

From equation ()105j) . we see that a is implicitly depends on r and consequently, from this 
equation it is not possible to predict the nature (decreasing or increasing) of dependence of 
a on r. But plotting a against r for the appropriate set of values of the parameters involved 
in the system, we find that a is slowly varying function of time. By the phrase “ appropriate 
set of values of the parameters”, we mean that those values of the parameters of the system 
for which the condition for existence of alternative solitary wave solution of the combined 
MKdV-KdV-ZK equation holds good, i.e., for those values of the parameters of the system 
for which L > 0. Taking ao = 0.5 (arbitrary) and the values of the parameters as mentioned 
in the figure [2l we plot a against r in Figj2j This figure clearly shows that the amplitude 
(a) decays slowly with time (r) and consequently, the amplitude of the alternative solitary 
wave solution of the combined MKdV-KdV-ZK equation is a slowly varying function of time 
when the effect of Landau damping is considered. 


V. CONCLUSIONS 

A macroscopic evolution equation corresponding to the combined MKdV-KdV-ZK equa¬ 
tion has been derived to include the effect of Landau damping. This macroscopic evolution 
equation admits the same alternative solitary wave solution of the combined MKdV-KdV- 
ZK equation except the fact that the amplitude of the solitary wave solution of the combined 
MKdV-KdV-ZK like macroscopic equation is a slowly varying function of time. The multiple 
time scale method of Ott and Sudan [3] has been generalized here to solve the said evolution 
equation. In small amplitude limit, we have observed the following result. 

Result:: Due to inclusion of the effect of Landau damping, the amplitude of the alternative 
solitary wave solution having profile different from seek 1 2 or seek of the macroscopic 
evolution equation decays slowly with time. 


[1] A. Bandyopadhyay and K. P. Das, Phys. Plasmas 9, 465 (2002). 

[2] J. Das, A. Bandyopadhyay, and K. P. Das, Phys. Plasmas 14, 092304 (2007). 


21 









[3] E. Ott and R. N. Sudan, Phys. Fluids 12, 2388 (1969). 

[4] G. W. Hammett and F. W. Perkins, Phys. Rev. Lett. 64, 3019 (1990). 

[5] J. F. McKenzie, F. Verheest, T. B. Doyle, and M. A. Hellberg, Phys. Plasmas 11, 1762 (2004). 

[6] T. S. Gill, H. Kaur, and N. S. Saini, Phys. Plasmas 10, 3927 (2003). 

[7] J.-K. Xue, Phys. Rev. E 69, 016403 (2004). 

[8] V. I. Karpman, J. P. Lynov, P. Michelsen, H. L. Pcseli, and J. J. Rasmussen, Phys. Fluids 
23, 1782 (1980). 

[9] R. A. Cairns, A. A. Mamun, R. Bingham, and P. K. Shukla, Physica Scripta T63, 80 (1995). 

[10] A. A. Mamun and R. A. Cairns, J. Plasma Phys. 56, 175 (2000a). 

[11] Y. Nejoh, Phys. Plasmas 5, 2830 (1992). 


22 







